Atomic Sole filter: multi-resonant photoemission via periodic poling of atom-cavity 

coupling constant 
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This paper describes a novel atom-cavity interaction induced by periodically poled atom-cavity 
coupling constant which leads to multiple narrow photoemission bands for an initially inverted 
two-level atom under the strong coupling condition. The emission bandpass narrowing has a close 
analogy with the folded Sole filter in the context of quasi-phase matching by periodic poling. We 
present a closed form solution of the emission probability at the end of interaction and deduce 
the multiple phase matching condition for this system which is programmable by the interaction 
time. The Bloch sphere analysis provides a clear understanding of the underlying atomic dynamics 
associated with the multiple resonances in the semiclassical limit. Furthermore, we show that this 
interaction can be applied to generation of the nonclassical field with sub-Poisson photon statistics. 

PACS numbers: 42.50.-p, 42.50.Pq, 42.50.Gy 
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I. INTRODUCTION 

Periodic poling of nonlinear medium [l], 0, 0, 3 is a 
widely used technique to achieve quasi-phase matching in 
nonlinear optics. The phase mismatch otherwise present 
in the medium is actively compensated by periodically 
reversing the phase slip so that the nonlinear response of 
the medium can be maintained over a large distance. The 
quasi-phase matching by periodic poling also finds its use 
in a polarization interference filter, namely the folded 
Sole filter 0,0]. It is constructed by periodically folding 
the fast axes of waveplates which are stacked between 
two orthogonal polarizers. A vertically polarized input 
is gradually brought into a horizontally polarized output 
for the quasi-phase matched angle of the fast axes. Since 
the intensity of the output is strongly dependent on the 
phase retardation, and in turn on the wavelength of the 
input field, the transmittance of the system is sharply 
peaked at target wavelengths. The passband of the filter 
gets very narrow as the number of stacks increases. 

In this paper we describe an atomic Sole filter, an 
atom-cavity system where a narrowing of photoemission 
band occurs via a novel atom-cavity interaction under 
the strong coupling condition. By periodically poling the 
atom-cavity coupling constant, one can actively guide the 
initially inverted atom into its ground state, thereby de- 
liver a photon into the cavity mode perfectly at certain 
atom-cavity detunings satisfying a quasi-phase matching 
condition. As the number of poling increases the band- 
width for the photoemission gets narrower. We present a 
closed form solution of the photoemission probability for 
an arbitrary number of poling. The quasi-phase match- 
ing condition, which is programmable by the atom-cavity 



interaction time, is deduced from the geometric consider- 
ation of the Bloch sphere. We explain the filtering mech- 
anism by examining the trajectories of the atomic state 
on the Bloch sphere. 

The atom traveling across a high-order transverse 
mode (TEMjvo) or along the cavity axis in a high-Q 
Fabry-Perot cavity may well be suited for the realization 
of the present idea. In particular, we show that the nar- 
rowing of the photoemission probability in the parameter 
space has a direct application to the generation of micro- 
maser/microlaser field 0, H| with enhanced sub-Poisson 
photon statistics [§, [l(| ■ 

This paper is organized as follows. In Sec. [Til a model 
Hamiltonian and a periodic atom-cavity coupling con- 
stant are introduced. Time domain analysis in Sec. IIIII 
provides the emission probability in a closed form by us- 
ing the unitary evolution matrix. In Sec. IIV1 the reg- 
ular structure observed in the emission probability and 
the phase matching conditions therein are interpreted in 
terms of the Bloch sphere. The analogy between the 
conventional Sole filter and our system is given in Sec. 
IVl Possible applications to nonclassical light generation 
and a frequency domain analysis are presented in Sec lVIl 
We conclude in Sec. I VIII 



II. MODEL 

We consider a two-level atom interacting with an elec- 
tromagnetic field mode of a cavity via electric dipole in- 
teraction. Under the rotating wave approximation and in 
the semiclassical regime of a large mean photon number 
n(^> 1) in the cavity mode, the interaction Hamiltonian 
is given by 
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FIG. 1: Periodically poled atom-cavity coupling constant 
shown as a function of time. 



where g is the atom-cavity coupling constant, a± are the 
Pauli pseudo-spin operators and A is the cavity-atom 
detuning. In our model system, the relative phase of the 
atom-cavity coupling constant is sequentially changing 
by 7r after each time interval of r, i.e., 



m = (-1) 



m+l 



.90 



for (to — l)r < t < mr , (2) 



where m=l, 2, 3, . . . , N (see Fig. [JJ. We call such an 
interaction the A^th-order multipolar interaction. The 
equation of motion is then 



c e (t) = -i^fng(t)t 



,(*) 



,(*) 



(3) 



(4) 



where c e (t) and c g (t) are the probability amplitude of 
the excited and ground state, respectively, in the inter- 
action picture. Our interest is to find the photoemission 
probability of an initially inverted atom after the entire 
interaction time of AY. 

Although the equations are written in the form of semi- 
classical Bloch equation, this model is also applicable to 
the vacuum Rabi oscillation in the restricted atom-field 
dressed basis of |e,0) and \g, 1) with sjn removed and 
2<7o understood as the vacuum Rabi frequency. Then the 
time evolution of the atomic system is exactly of the same 
form for both quantum and semiclassical fields if we trace 
out the photonic variables in the quantum calculation. 



III. EMISSION PROBABILITY 

If we denote the atomic wave function in the inter- 
action picture as a column vector, the evolution of the 



atomic state is described by a unitary matrix as 
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where ci" 1 -* and Cg ^ are the excited and ground state 
probability amplitude at t — mr, respectively. The ele- 
ments of C/( m ) are given by 



(U(m)\ = { _ iri H sin ^(m-l/2)S 
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where i] = 2y/ngor, S = At, and (f> = \J "g 2 + S 2 , the Rabi 
precession angle. As mentioned in Sec. [Til this matrix 
is same for the vacuum Rabi oscillation and the semi- 
classical one. Note also that we should keep the index 
to indicating the sequential order in order to incorpo- 
rate the correct initial conditions for each stage since 
[fr/(ti),fl>(* a )] ^Oingeneral. 

We are now interested in the probability by which an 
initially inverted atom emits a photon into the cavity 
mode. The probability, called photoemission probability 

P em , is nothing but the ground state probability Ic^'l 2 
of the atom after the entire interaction time Nr. For 
calculating P em it is convenient to define a unit cell op- 
eration which is comprised of two intervals with ±go as 



T (k) = u (2k)jj(2k-l) 
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where its elements are 
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For even N the final state is then simply 
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TABLE I: Emission probabilities after the entire interaction time Nt 
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Note that each matrix element of is independent 
of the sequential order k except the extra phase factors 
e ±2ife<5_ -yy e can s jj 0W ^hat those phase factors do not 
affect the absolute value of the probability amplitude (see 
Appendix A) so that we may eliminate them from 
and define the resulting matrix T as 
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This elimination of fc dependence reduces the complexity 
of the calculation much. From the Chebyshev identity 
[ill ] the product of an unitary unimodular matrix is given 
as 
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where cos£ = ^(A + D). Hence, the final results for an 
initial state 
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FIG. 2: Density plot of emission probability in S-r/ plane. 
From (a) to (h) N increases from 1 to 8. The dashed lines 
(called branch lines) indicate the maximum visibility of Rabi 
oscillation and designated by the index p in the text. We 
can also recognize the semicircular patterns (called branch 
circles), which are designated by q. Quasi-phase matching is 
accomplished at the intersections of the branch lines and the 
branch circles. 



Using the Chebyshev polynomial in the trigonometric 
form [i~2| one can expand 
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to have the emission probability as a function of Rabi 
frequency (77) and atom-cavity detuning (6). First few of 
them are explicitly listed in Tab[T]and plotted in Fig. [2j 
The result for N = 1 is the familiar Rabi oscillation 
(Fig. [2Ja))- The visibility of the oscillation is maximum 
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on resonance (S — 0). For N > 1 the additional regular 
patterns emerge in the (77, 6) map as in Fig. [Ub)-(h). In 
particular, the 100% certain emission of a photon is pos- 
sible for off-resonant conditions. We notice full visibility 
along a specific straight line, i.e. at a fixed ratio of r\ 
to 5. Furthermore, the fringe pattern gets sharper for 
the smaller r)/5 ratio. In the next sections we investigate 
what physical processes are involved for such patterns 
and the sharpness of the fringes. 



IV. BLOCH SPHERE ANALYSIS 

The coherent evolution of the internal state of a two- 
level atom is easily visualized as the rotational motion of 
a state vector around a torque vector on the Bloch sphere 
in the semiclassical limit [13j ■ The equation of motion is 

R = Rxft (18) 

where R = [§ (c e c* g + c g c* e ) , ■§ (c e c* g - c g c* e ) , |c e | 2 - |c 9 | 2 ] 
the state vector, and fl = [2y/ng(t), 0, A] the torque 
vector. Initially the state vector R is pointing the north 
pole. When the interaction is turned on, R precesses 
around the torque vector Q. Flipping the polarity of the 
coupling constant is equivalent to fold the torque vector 
with respect to z-axis in the x-z plane (as in and 
O in Fig. [4]). 

Suppose the cavity is tuned on resonance, correspond- 
ing to the line of 8 = in Fig. [2] Since $1 is pointing the 
exactly opposite directions periodically, the state vector 
just retraces a part of the great circle as in Fig. (3jb). 
Thus for even N the interaction ends up with no emis- 
sion of a photon. We see the emission probabilities are 
completely quenched along 5 = in Fig. [2] for even N. 
For odd N, on the other hand, the interaction is effec- 
tively unipolar since after (N— l)r the state vector comes 
back to its initial position. Indeed along S = the emis- 
sion probabilities are the same for all odd N as shown in 
Fig.H 

For off-resonance, however, the multistep flipping of 
the torque vector brings about various trajectories on 
the Bloch sphere. Some of such nontrivial motions of the 
state vector lead to complete emission of a photon even 
at off-resonance unlike the conventional unipolar interac- 
tion. In the S-rj plane of Fig. [2] these occur only along 
some specific ratios of i]/S, i.e. at specific angles of the 
torque vector. Some of such trajectories of the state vec- 
tor are shown in Figs.[3£c) and[3Jd). The trajectory with 
the smallest rj/8 [Fig.^d)] is noticeable among them in 
that it steadily goes down to the ground state, utiliz- 
ing all N steps of interaction in increasing the emission 
probability. In the other trajectories, on the other hand, 
the interactions in some intervals are contributed to ex- 
cite the atom back although the state vector eventually 
arrives at the ground state in the end. 

The requirement for the complete emission, or namely 
the quasi-phase matching condition, can be found via 




FIG. 3: (a) Emission probability near = 7T (q = 1) semicir- 
cular branch for N — 5. The red crosses indicate the points 
for the complete emission of a photon. The corresponding 
trajectories of the Bloch vector are shown for (b) p — 3, (c) 
p — 2 and (d) p = 1. The red curves represent the evolution 
under +go while the blue curves are under —go as in the same 
color convention of Fig. [T] 




FIG. 4: Evolution of the Bloch vector in x — z plane assuming 
(j> — 7r. The filled circles are the location of the Bloch vector 
after each r and fl± are torque vectors associated with ±go, 
respectively. 

a simple geometric argument. Assume the accumulated 
net precession angle <f> in Eq. ([6|) is odd multiples of tt for 
each interval. This condition constitutes the semicircles, 
called branch circles, in the map of P em in Fig. [2] Then 
the atomic state after each interval reduces onto the x-z 
plane of Bloch sphere as in Fig. |H where the position 
of the state vector is numbered for a few sequence of 
interaction. 

In (x, z) coordinate a reflection about the torque vector 
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is represented by 
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for (m — l)r < t < mr with 9 being the angle between z- 
axis and the torque vector. Then the unit cell operation 
of N = 2 is represented by 



cos 49 -sin 40 
sin 49 cos 49 



(20) 



and the generalization to an arbitrary even N is straight- 
forwardly given as 



cos2N9 -sm2N0 
sin2N9 cos2N9 



(21) 



Applying one more reflection to Eq. (j2"Tj) gives the odd N 
result which turns out to be 



-cos2A6> sm2N9 
sin 2N9 cos 2N9 



(22) 



Therefore the position of the state in (x, z) coordinate 
changes as 



(23) 



Thus the quasi-phase matching occurs for cos2N9 = — 1, 
that is, 
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where p = 1, 2, • • • , N. These particular angles constitute 
the straight lines, called branch lines, in Fig.[5J Since r/ — 
<5tan6> and <f> = y/rj 1 + 8 2 = (2q— l)ir, where q = 1, 2, • • • , 
the optimal detuning for complete photocmission is given 
by 



Sopt/n = ±(2q- 1) cos 



N 



(25) 



where p and q index the branch lines and the branch 
circles, respectively. 



V. SOLC FILTER ANALOGY 

It is noted in Fig. [2] that the widths of fringes along 
different branch lines are significantly different from each 
other. The fringes are sharper along a branch line with 
smaller \9\. As mentioned in the previous section, the 
state-vector trajectories associated with different p's are 



qualitatively different from each other. In particular, pre- 
cise multiple rotations are sequentially involved for the 
trajectory of Fig. [3f d), and thus a slight slip from the op- 
timal condition will lead to a significant reduction in the 
photoemission probability. That is why the fringe along 
the branch line of p = 1 (or equivalently p — N) is the 
sharpest. Furthermore, the passband of 6 gets narrower 
as N gets large. 

This behavior is analogous to what happens in the op- 
tical Sole filter, where the waveplates with their fast axes 
inclined by ±9 about the input vertical polarization are 
periodically stacked and the phase retardation angle (f> 
imposed by each waveplate is set to n (half wave plate) 
for a target wavelength. Here the angles 9 and <j> in the 
optical Sole filter work as nothing but the inclined angle 
of the torque vector and the net precession angle, re- 
spectively, in our atom-cavity case while the two kinds of 
polarization correspond to the two energy states of the 
atom and the Jones matrix of the filter corresponds to 
the matrix T. 

In fact, one-to-one correspondence can be explicitly 
shown if we restrict ourselves to the absolute value of 
the probability amplitude c ej g(7Vr). If we perform the 
coordinate transformation 



tan 1 {r]/8) 



5 2 



(26) 



we can obtain the spectrum of filtering (actually pass- 
band of photoemission) of our atom-cavity system as a 
function of <f> for given azimuth 9 as in the optical Sole 
filter. We confirm two aspects of filtering functions in 
Fig. El Firstly the passband gets narrower for larger N 
and secondly for given N the lowest p (or equivalently 
the largest p) gives the sharpest passband of photoemis- 
sion. Due to this analogy our atom-cavity system can be 
viewed as an atomic Sole filter. 

However, we also recognize the difference between the 
two systems when we take the phase of the probabil- 
ity amplitude into account such as the atomic polariza- 
tion, for which the optical Sole filter has no correspond- 
ing counterpart. Note then we have to use the phase 
sensitive matrix T^ k > instead of T to retrieve the phase 
information in our atom-cavity system. 



VI. DISCUSSION 

A. Sub-Poisson photon statistics 

The squeezed lineshape of P em in S — i] plane, brought 
about by the multi-polar interactions, is useful for gener- 
ation of nonclassical field state with sub-Poissonian pho- 
ton statistics in the cavity-QED micromascr/microlaser. 
While the mean photon number of the conventional laser 
is stabilized by the cavity damping only [Io| . the active 
role of coherent medium is also very important in the 
cavity-QED lasers. Specifically, the negative slope of P em 
as a function of the mean photon number n reduces the 
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FIG. 5: (a) For given phase matching order p — 1 (or t9 = 
tt/2N), the passband of the atomic Sole filter, our atom-cavity 
system, is shown as a function of <f>(r], 8). (b) For given number 
of poling N = 8, the passband for different phase matching 
order p are shown. 



recovery time of the system when the intracavity pho- 
ton number is momentarily deviated from its stationary 
point [id ], and thereby can reduce the photon number 
fluctuation below the shot-noise limit, resulting in a sub- 
Poissonian photon statistics. On the contrary, the slope 
of P em is always positive in the conventional laser. As a 
measure of the photon number stabilization, we use the 
Mandel Q factor, which can be written as 



dPem/dn 



dP em /dn 



(27) 



where D cav is given by the cavity decay rate divided by 
the atomic flux [id] . The negative value of Q indicates a 
sub-Poissonian photon statistics. 

We plot the Mandel Q in Fig. [H] for various N in n — 5 
plane with 770 = 2goT fixed at 0.24 and -D cav =0.0038, both 
of which are typical experimental parameters of Ref. [Io| . 
Due to the squeezed structure of P em for large N's, we 
expect the more negative slope of P em for the higher- 
order multipolar interaction. However, the direction of 
squeezing is along the branch lines, not along n, except 
for the first branch circle, along which the squeezing is 
also pronounced as shown in Fig. O Therefore, the nar- 
rowing effect along n is only good for small detuning of 
5/tt ~ ±1 for large N. For large detunings the direction 
of squeezing is roughly along S for the first and the last 
branch lines (p — 1,N). We find the first/last branch 
line solution with <f> = 7T (q = 1) provides the most nega- 




FIG. 6: Mandel Q factor for 770=0. 24 and D C av=0.0038 as N 
is varied from 1 (a) to 8 (h). Only the negative part is shown 
for simplicity. 



tive slope of P em along n for larger N, and thus the most 
negativity of Q(~ —1) or the most nonclassicality of the 
field as shown in Fig. [S] 

In an experiment employing an atomic beam with a 
velocity spread, the slope of P em would be reduced due to 
velocity averaging. However, the peak giving rise to the 
most negative slope (p — {l,N},q = 1) is least affected 
by this averaging and thus its Mandel Q is expected to 
be somewhat robust to such statistical averaging. 



B. Fourier analysis and large iV limit 

The coupling constant in Eq. ([2]) is a truncated peri- 
odic square-pulse train and can be expressed by a Fourier 
series on the time interval [0, Nt] as 



i{t) = £ 



U/Nt 



(28) 
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where 



J(0 



At 



(29) 



Likewise the probability amplitudes can be written in the 
similar form 



(l) -Vhrlt/Nr 



(30) 
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The differential equations Eq. © and Eq. ^ are trans 
formed into algebraic equations of Ce} g as 



C (i) = NV V- 



-k±s) 



k— — oo 



2.9o 



_ r (fc) 



(31) 



respectively, where we index the detuning with an integer 
s as A = 2tts/Nt. Initially only Ce is nonzero. As 
the interaction proceeds the Fourier elements G"' of the 
coupling constant generates the ground state probability 
amplitude. 

For small r\ < 1, so called the weak coupling limit, 
the Fourier coefficients are coupled approximately in the 
first order of 77. Since the precession angle is very small, 
Ce^ almost stays in its initial value and d is dominant 
over the summation in Eq. [)3ip. Then Cg , which is 
the dominant frequency components of the ground state 
amplitude, resonantly grows for the frequencies at which 
G^ is significant. Therefore, the ^-dependence of P em 
almost follows the spectral lineshape of git) in the first 
order approximation. Further occupation of Cg is ac- 
complished for the strong coupling condition (77 > 1) and 
the rich structure of P em is then revealed. 

In the limit of large N, the coupling constant is well ap- 
proximated by the sinusoidally modulated function which 
in turn leads to the bichromatic atom-field interaction. If 
the modulation frequency (tt/t here) is much larger than 
the Rabi frequency, the system is well described by the 
shifted resonance [14]]. In fact, for rj <C 1 the emission 
probability is peaked at S ~ ±7r, which are shifted reso- 
nance frequencies. This claim becomes more correct as 
TV increases as we see at the very bottom of each map in 
Fig. [2] because the coupling constant is more like a sinu- 
soidal function. For 77 > 1 it is known that the harmonic 
resonances come in to play so that multiple branches of 
peaks in the emission probability emerges as in Fig.O Al- 
though the bichromatic interaction has been extensively 
studied [3, [H, EE E3 j tne systematic analysis for the 



finite train of flipping Rabi frequency as given here never 
appeared before. 



VII. CONCLUSION 

We have studied the atom-cavity system whose cou- 
pling constant is periodically poled N times in time. We 
have shown that the photoemission probability of an ini- 
tially inverted two-level atom at the end of the inter- 
action exhibits narrowed passbands determined by Sole 
filter-type quasi-phase matching condition. The closed 
form expression for the emission probability as a func- 
tion of Rabi frequency and detuning is provided for an 
arbitrary number of poling. From the Bloch sphere inter- 
pretation we can verify the underlying physical processes 
of the quasi-phase matching in the emission probabil- 
ity. The narrowing of the photoemission probability in 
the parameter space can be applied to the generation of 
micromaser/microlaser field with enhanced sub-Poisson 
photon statistics. 
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APPENDIX A: PHASE FACTORS IN THE UNIT 
CELL OPERATION 

The transfer matrix is a unitary unimodular ma- 
trix. Along with the particular form of its phase factors 
it can be written in its most general form as 



T (k) = 



ae l6 
_ b * e -i(2k- 



bet (2k-l)8 



(Al) 



where a and b are arbitrary complex numbers satisfying 
the unimodular relation \a\ 2 + \b\ 2 = 1. The product of 
such matrices is of the form 



N 

n r(fe) 



*:=i 



fia,b)e mf> g(a,b)e im 
-g*(a,b)e- iNe /*(o, b) e - me 



(A2) 



where the functions /(a, b) and g(a, b) are some polyno- 
mial functions of a and b. The emission probability cor- 
responds to \g(a, b)\ 2 , in which the common phase factors 
do not appear. 
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